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ABSTRACT 

We generalize the mean field magnetic dynamo to include local evolution of the 
mean vorticity in addition to the mean magnetic field. The coupled equations exhibit 
a general mean field dynamo instability that enables the transfer of turbulent energy 
to the magnetic field and vorticity on larger scales. The growth of the vorticity and 
magnetic field both require helical turbulence which can be supplied by an underlying 
global rotation. The dynamo coefficients are derived including the backreaction from 
the mean magnetic field to lowest order. We find that a mean vorticity field can actually 
seed exponential growth of mean magnetic field from only a small scale seed magnetic 
field, without the need for a seed mean magnetic field. The equations decouple when 
the fiuctuating velocity and magnetic field cross-correlations vanish, resulting in the 
separate mean field magnetic and mean field vorticity dynamo equations. 



Subject Headings: magnetic fields: MHD; instabilities; accretion disks; sun: magnetic fields; galaxies: 
jets. 
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1. Introduction 

The mean-field magnetic dynamo (MFMD) 
has been a standard framework for understand- 
ing the origin of large scale magnetic fields (B- 
fields) in planets stars and galaxies (Moffatt 
|97|, [Parker 1979| , [Ruzmaikin et al. 1988| , |Beck 



t al. 1996). In this mechanism, a large scale seed 



field grows exponentially at the expense of small 
scale turbulent energy. This can be distinguished 
from fast dynamo turbulent amplification of B- 
fields (e.g. Parker 1979| ), which occurs on the 
scale of the input turbulence. Many systems 



such as planetary and solar spots (Peter 1996) 



and accretion and galactic disks ( Abramowicz 



dei LOC for large ocalo vortcic otruotufGO. Jcta may 



19(^ — Kitchatinov et al. 1994a) also show ovi 



tially zero, in contrast to the decoupled theory. 

We include the mean B-field backreaction on 
the velocity flows to lowest order, thus deriving 
corrections to the kinematic magnetic and vor- 
ticity dynamo coefficients. However, we do not 
compute these corrections to higher order, nor 
do we purport to offer a fully developed non- 
linear theory of MHD turbulence. We make ex- 
plicit assumptions/approximations, but empha- 
size that none of these are beyond those used 
in standard mean-field dynamo treatments. (In 
many treatments, the same approximations are 
imposed without explanation.) We are aware 
of the controversies of the non-linear backreac- 
tion (e.g. Piddington 1981 , Kulsrud &: AndeF| 
son 19821, [Field fc Blackman, 199^, [Cattaneol 



als' J bu iulalud. In Lliu abtjuncu uf B fiuld i ^j, t Jrir 



1994, Pouquet et al. 1976, Blackman 1996, Bran 



me an vorticity equation is similar to the mean 



magnetic induction equation (Kitchatinov et al. 



denburg 1996| , pubramanian 199^ ) , but inasmuch 
as the "a^"(e.gr iM^ffatt 1978D or "a - 17" (e.g. 



1994a]) and can allow mean field vorticity dy- 
namo (MFVD) growth. When the mutual cou- 
pling of vorticity and B-field is considered, the 
mean quantity whose time evolution is of inter- 
est has 6 components: 3 components each for the 
mean B-field and mean vorticity. Growth arises 
in this coupled mean field dynamo (CMFD). 

We first derive the equations for mean and 
fluctuating quantities and show that the evolu- 
tion equations for the mean vorticity and B-field 
are non-trivial only when the turbulence is at 
least weakly inhomogeneous and anisotropic. By 
expanding the turbulent quantities to first order 
in the time-varying mean velocity and mean B- 
field, we derive the initial mutual dynamo growth 
for several sets of initial conditions. When cross- 
correlations between turbulent velocities and tur- 
bulent B-fields vanish, mutual field growth is de- 
scribed by decoupled MFVD and MFMD equa- 
tions. However, when the cross correlations do 
not vanish, the time evolution of the mean fields 
is coupled. As a result, we find that a seed mean 
vorticity and small scale B-field can allow growth 
of the mean B-field even when the latter is ini- 



Parker 1979| ) dynamos are at least a framework 



for the study of B-field generation and compar- 
ison with observations, so too is the more gen- 
eralized treatment of the vorticity-magnetic field 
dynamo developed herein. 

2. Derivation of Coupled Mean-Field Equa- 
tions 

Consider the Navier-Stokes equation for in- 
compressible flows in the presence of B-fields, 



= V X (V X v) - Vpeff + i^V^v b • Vb 

+F(x,t) + V(/>, 
(1) 

where is a constant viscosity, F is a forcing 
function, p^^jj = p + 6^/2 -|- w^/2, with p = P/p 
and b = B/(47r/9)^/^ is the constant density nor- 
malized B-field. The V(/) includes potential forces 
such as gravity. The vorticity equation. 



dtu; = Vx(vxu;)-hi/V^a;+Vx(b-Vb)-hVxF(x,t), 

(2) 
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where a; = V x v, is obtained by taking the curl 
of dll). The induction equation for b is 



5tb = V X (v X b) + z^AfV^b, 



(3) 



where vm is a constant magnetic viscosity. 

Next, decompose u:, v, and b into mean (in- 
dicated by an overbar or ( )) and fluctuating 
(indicated by prime) components v = vj- + v', 
u} = u^T + and b = b + b', respectively with 
= V + Vc- Here, Vc is a constant global ve- 
locity whose curl is a>c) and v and a; are the time 
dependent mean components whose spatial scales 
of variation are smaller spatial than that of a>c- 
The u3c supplies helicity to the turbulence and 
ensures that any growth of oi represents a trans- 
fer of angular momentum from the turbulence. 
(Alternatively, the helicity could be supplied by 
F(x,t).) 

We assume that derivatives with respect to 
x_or t obey dt,x{XiXj) = {dt ,x{XiXj)) an d 
{XiX'j) = (Reynolds relations (|Radler 1980| )), 
where Xi = Xi + X[ are components of vec- 
tor functions of x and t. For statistical ensem- 
ble means, these hold when the correlation time 
scales are short relative to the variation times of 
mean quantities. For the spatial mean, defined 
by (Xi(x,t)) = |Cr3 4"+/'Xi(s,t)d3s, the rela- 
tions hold when the averaging is taken over a 
large enough scale, such that / <C |C| ^ where 
L ~ 6/V6, and I ~ h' /Vh' ~ v'/Vv'. 

Subtracting the mean of (|l|) from itself, and 
assuming V(/>' = 0, gives 



d^ y' = (v' . Vv^) - . Vv^ - V . Vv^ - . Vv- 



(b' • Vb') + b' • Vb + b • Vb' + F'(x, t) + z^VV, 

(4) 

where dt = dt + Vc • V . The average of (|2|) is 
dttD = V X (v' X w') + z^V^WT + V X (b' • Vh') + 

i^c ■ Vv + U3 ■ VVc + U}c ■ VVc, 

(5) 



where we have neglected terms second order in 
time-varying mean quantities. Subtracting (^) 
from (^) gives 

dfu^' = bJ ■ Vvy — V • u;' + (It • Vv' — v' • Vw' 
-v' • VwT + ^' ■ Vv' - V X (v' X J) + i/V^w'-F 

V X (b' • Vb) + V X (b • Vb') + V X (b' • Vb')- 

V X (b' • Vb') V X F'(x,t). 

(6) 

Similarly, the equation for the mean B-field, de- 
rived by averaging (^ is 

dtb = V X (v' X b') + b • Vvc + z^AfV^b. (7) 

Subtracting (|^) from (|3|) yields the equation for 
the fluctuating B-field 

dtb' = b' • VvT - V • Vb' b • Vv' - v' • Vb+ 
b' • Vv' - v' • Vb' - V X (v' X b') + z^A^V^b'. 

(8) 

Growth of b or requires the Vx terms in 
(|7|) and (||) to be non- vanishing. However, when 
the turbulence is strictly homogeneous in v' and 
b', these quantities vanish. For purely isotropic 
turbulence, the term (v'xb') in ^ vanishes since 
it is the average of a vector, while V x (v' x u;') oc 
V X V(f;'2) and V x (b' • Vb') oc V x V(6'2) vanish 
from Reynolds rules and incompressibility. Thus, 
anisotropy and inhomogeneity must be present 
for nontrivial time evolution of mean fields. 

We expand the turbulent quantities on the 
right of Eqs. (^) and (^) to linear order in hoih 
b and v using the equations for the fluctuating 



fields thus generalizing the approach of (Field & 



Blackman, 1996 ) where only the mean B-field 



was used. To find the lowest order terms, we 
assume weakly anisotropic inhomogeneous tur- 
bulence: Terms linear in the time-varying mean 
quantities contribute, but their averaged 0*^ or- 
der coefficients, are taken to be isotropic and ho- 
mogeneous (still allowing for reflection asymme- 
try) . Iterating the equations using the formal so- 
lutions for the turbulent fields b'(t) = b'(t = 0)-|- 
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/ dt>h'dt' and = a;'(t = 0) + / dfJdt', and 

using times appropriately chosen such that the 
correlations (v'(t) x a;'(0)) = (v'(t) x b'(0)) ~ 0, 
we obtain to first order in mean quantities 



(v'xa;')« = (v'W(t)x/o*d,,u;'WdO+ ... 

(/o*d,,v'«rft'x.^'(0)(t)), 

wi th similar expressions for (b^ ■ VbO(^) and 



where V ■ Vc — 0, but V • v, V • v' 7^ 0. The aver- 
aged terms in (^l]) make no contribution when in- 
serted into ( |lO| ) according to the Reynolds rules. 
We only consider terms explicitly linear in time- 
varying mean quantities, and iterate to first or- 
der in the correlation time; higher order terms 
would come from successive iterations in the 
equations for the fluctuating variables. This first 
order smoothing approximation (FOSA) ( [Krause 



(v ix b')(^. The calculation of these averages & Riidler 198O0 is "justified" when the correla 



requires Eqs. 
gr ands. Using 
thd 



via 



(^), and (§) for the time inte- 
also requires an expression for 



pressure, which arises in (^1 and in (v^ x b^) 



the terms 



(a;'(0)(t)x / Vp'^^Ut') and (b'(o)(t)x / Vp'^^^dt') 
JO Jo 

(10) 

Using isotropy, homogeneity, and Reynolds 
rules, we now show that remarkably, terms of the 
form (|l^) vanish in the derivation of the mean 
field equations to the order considered. Con- 
sider the pressure in the isoentropic energy equa- 
tion (hereafter neglecting microphysical fluid and 
magnetic viscosities) DtP = (itP-|-(v'-|-v)-VP = 
-7PV • V, where 7 is the adiabatic index ([Batch 



elo j: 1967 , Boyd and Sanderson 1969| ). We have 
assumed incompressibility (V • v ~ and con- 
stant p) above, but here we must allow V ■ v 7^ 
for the last term in the energy equation: Even 
though V • V is small, the sound speed is large 
(i.e. 7 is large, see Patchelor 1967 ) such that the 
magnitude of DP/Dt cannot be ignored. This is 
why the energy equation is normally not useful 
in incompressible approaches. But here it will be 
symmetry properties of the energy equation that 
account for the absence of a pressure contribu- 
tion when used in (^) and the seemingly offend- 
ing term will not contribute. Taking p' = P' / p 
we obtain 

dtp' = (v' • Vp') - v' • Vp - V • Vp' - v' • Vp'+ 
7(p'V • v') — 7^)' V • V — 7j5V • v' — 7^ V • v', 

(11) 



tion time of the turbulence is less than the eddy 
turnover time ( Ruzmaikin et al. 1988| , Blackman 



1995), as is possibly the case in galaxies ( [Ruz- 



maikin et al. 1988|) , or if higher order correlations 



are otherwise reduced. Integrating (pT[), and us- 
ing the assumptions discussed above. 



V/^) = v/^)(0) -V j (v'(°)-Vp + v-V/°) (12) 
+^/0)V-v + 7pV-v'(°))dt'. 

The pressure dependent contribution to (v'^'^) x 
u}'^'^'^) can then be written 

rt 

X / 4Vp'dt')k = eijkeimsX 



'dt' f dt"{{djvidip'^'^+vidjdip"^'^^ + 
Jo 

djv'i^^dip + v'l^^djdip + ^p'^°^djdivi+ 
jdjp'^^^ divi + 'ypdjdiv'i^'^^ -\- jdjpdiv'^^^^) 

= i f'dt' f dt"{J^'\t').J^^){t"))d,p{t), 
6 Jo Jo 

(13) 

where we assume mean fields vary on time scales 
longer than the fluctuating fields, and \/p'^^\t = 
0) is uncorrelated with b'(i) or u:'{t). The only 
surviving term arises from the third term on the 
right hand side of (^3|) and the vanishing of the 
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remaining terms follows from careful application 
of isotropy {i.e. rank 2 and rank 3 averaged ten- 
sors of fluctuating 0*'^ order quantities are pro- 
portional to 5ij and €ijk respectively) and homo- 
geneity {i.e. = 0) of the 0*^ order 
turbulence, without invoking V • v' = 0. The 
surviving term in (|l3|) vanishes when put inside 
the curl in (^). A similar analysis holds for the 
second term in ( |lO|) when put into (0); thus, the 
pressure does not contribute to (|5|) or (|7|). 

Approximating time integrals by factors of the 
correlation time Tc ( Ruzmaikin et al. 1988| ), and 
freely employing Reynolds rules and incompress- 
ibility we obtain 



(3/t,)V X (v' X a;^^) = ^c^f ^)(V x u>) + 

(e.,..f9.6f)(Vxb)-(.f)6:(°))V2(Vxb) 

(14) 

with similar expressions for Vx (b'-Vb')(i) and 
V x (v' x b')(^\ Upon substituting these into 
@ and (0), the curls can be pulled onto the u; 
and b from homogeneity of O*'^ order averages. 
Note that a>c and Vc contribute to the 0*'' order 
quantities, and thus do not show up explicitly in 
(13). After some simplification, we obtain 

dtti; = ao(V X w) + ai(V x b) + /^oV^tl; + /3i V^b 



-a2V2(V x b) + cZ>c • Vv + tl; • Vvc + ■ Vv, 



(15) 



and 



dth = a2(V X w) + a3(V x b) + /32V^b + b • Vvc 

(16) 

where the coefficients are: 



ao = (re/3)(a;'(o).v'(o)) 

ai = (re/3)(a;'(0) . V X b'W) 

02 = (2rc/3)(b'(0) • v'W) 

as = (2T,/3)(b'W • V X b'(0)) - ao 

/?o = (re/3)(v'(0) . v'(0) + b'(0) • b'(0)) 

/3l = (2r,/3)(a;'(0) -b'W) 

/32 = (r,/3)(v'(0) • v'(0) + 2b'(0) • b'(o)), 

3. Solutions and Discussion 



[17) 



The coefficients as and P2 in ( [Tq ) are modi- 
fied helicity and diffusion coefficients of the stan- 
dard magnetic dynamo theory. A similar, but 
not identical, set of helicity and diffusion coef- 
ficients enter ( [l5| ) for the mean vorticity. The 
cross-correlations 01,02,(^1, taken to vanish in 



Vainshtein & Kitchatinov 1983, need not vanish, 



as they preserve all MHD symmetries in ( [L5|) and 
([T^. We will now see how the finite helicity is 
required for both B-field ( [Moffatt 1978| ) and vor- 
ticity growth ( Kitchatinov et al. 1994a| ). 

To show the growth, we assume constant dy- 
namo coefficients and consider the simple case 
when the last 3 terms in ( p!5[ ) and the last term 
in ( [T^ ) are negligible. This assumption is valid in 
systems for which a>c is perpendicular to the gra- 
dient of the time varying mean quantities, and for 
which the former is relatively uniform or perpen- 
dicular to the dominant B-field component. (The 
influence of is always present however, in pro- 
viding reflection asymmetry to the 0*^ order tur- 
bulent quantities.) The simplified equations yield 
a generalized vorticity-magnetic a^ dynamo. 

The mean-field equations are then solved by 
assuming solutions of the form exp{at — ik ■ f) 
to obtain the linear algebraic equations aX = 
M.X , X = {Ox ,ujy,Ldz,bx,by,bz). The eigenval- 
ues a of the dynamical matrix M(A;) determine 
the time evolution of the fields at wavevector 
k. Different behavior can be inferred from the 
general solution X{f,t) = ^j^i^^Cj^j exp{ajt — 

ik ■ r). Here, Cj are depend on the initial val- 
ues Xj{t = 0), and are eigenvectors which are 
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linear combinations of ui and b. When cross- 
correlations between the turbulent velocities and 
B-fields are non-vanishing, the mean vorticity 
and B-fields are naturally coupled. The general 
6x6 matrix equation can be solved numerically 
for Xj and contains all decaying and growing 
modes. Because the Cj are combinations of the 
initial Xj(0), for growing modes, all components 
of b and ui will in general grow unless very spe- 
cific initial conditions X(t = 0) and values of k 
and (T{k) make certain Cj vanish. 

First we consider thin disk systems which have 
only in-plane mean field components (while al- 
lowing the turbulence to be 3-D), i.e., v = 
{vx,Vy,Q),h = {bx,by,0). The three possible 
poles are ao = — /JqA;^ and 2a± = — (/3o + ^2)^'^ i 



|fc|y4aia2 + iPo - l32Yk'^ + 'ia^k'^. which have 
growing modes {a > 0) when k'^ < 0102/ {/3oP2 — 
02) (if and only if aia2/(/3o/?2 — "2) > 0)- Thus, 
the growth of tOz and (bx, by) depends on the signs 
and magnitudes of ai,a2, Po, and 02- because we 
allow the small scale 

The fully 3-D case requires numerical calcu- 
lati on of (7j., but wo can obtain approximate bo 



ha vior by truncating ( tS ) and ( 3:6 ) to lowest or- 
der in mean field gradients. By taking the curl 
of ( p!6| ) and using j we see that terms with one 
higher derivative in the mean quantities are re- 
duced by a factor i/L. If 02 oc (v'(°) • b'(o)) 
is negligible initially, it will remain so ( phan- 
drc.n 1997, Kraichnan 19591) . Neglecting only 

the six roots are 
-P2k'^ ±a3\k\, 



the a2 term in ( |15| ) and ([Tf_ 

a = -Pok'^, -P2k'^, -Pok"^ ± ao|A:| 
which displays two growing modes for small k de- 
pending on the signs of ao and {Po, Ps > 0). 

Note that if cross-correlations between turbu- 
lent velocities and magnetic fields vanish, the 
equations for u) and b decouple into separate 
MFMD and MFVD equations of the form dttn = 
aV X m + /3V^m where m = or b, with their 
respective helicity and diffusion coefficients. The 
time evolution of the decoupled dynamos is de- 
termined by the eigenvalues a = —Pk'^, ±a|A;| — 



Pk"^. Thus, there is a growing solution (cr = 
\a\\k\ — pk"^ > 0) for small wavevectors \k\ < 
\a\/p. The most unstable mode is k* = |a|/(2/3) 
which grows at a rate a{k*) = ja^|/(4/3), like the 
usual dynamo ( [Moffatt 1978|) . 

To conclude: We have generalized the mean- 
field dynamo to include mean vorticity and have 
considered solutions of an vorticity-magnetic 
dynamo. An important implication of (|l^) and 
([1^ ) is that the mean B-field can grow even when 
its initial value is zero because a seed mean B- 
field can be generated by a coupling between the 
small scale B-field to the small scale velocity, and 
the latter to the mean vorticity. This feature 
should be preserved in more general treatments 
as it simply results from the coupling between 
( [T5|) and ([T^ and only requires non-vanishing 
cross correlations. Growth requires helical tur- 
bulence which can be supplied by an underly- 
ing global rotation. Astrophysical applications of 
more extensive treatments may include accretion 
disks, where vortices and B-fields can combine 
to concentrate emerging radiation into strong 
beams ( lAcosta-Pulido et al. 19"90| , ^ ramowic^ 
"T992D . Also, Yoshizawa & Yokoi (1993) showed 



that an equation similar to (|l^) can lead to 
generation of toroidal fields whose upward pres- 
sure drives a jet that subsequently may develop 
poloidal fields. The CMFD may provide the sus- 
taining feedback. 
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